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ABSTRACT 

We report on our analysis of the 1 Ms Chandra observation of the supernova remnant Cas A in 
order to localize, characterize and quantify its non-thermal X-ray emission. More specifically, we 
investigated whether the X-ray synchrotron emission from the inside of the remnant is from the 
outward shock, but projected toward the inner ring, or from the inner shell. We tackle this problem 
by employing a Lucy- Richardson deconvolution technique and measuring spectral indices in the 4.2-6 
keV band. 

We show that most of the continuum emission is coming from an inner ring that coincides with the 
location of the reverse shock. This inner ring includes filaments, whose X-ray emission has been found 
to be dominated by X-ray synchrotron emission. The X-ray emission from these filaments, both at 
the forward shock and from the inner ring, have relatively hard spectra with spectral index > —3.1. 
The regions emitting hard X-ray continuum contribute about 54% of the total X-ray emission in the 
4.2-6 keV. This is lower than suggested by extrapolating the hard X-ray spectrum as measured by 
BeppoSAX-PDS and INTEGRAL. This can be reconciled by assuming a gradual steepening of the 
spectrum toward higher energies. We argue that the X-ray synchrotron emission is mainly coming 
from the Western part of the reverse shock. The reverse shock in the West is almost at rest in our 
observation frame, corresponding to a relatively high reverse shock velocity of ~ 6000 kms -1 in the 
frame of the freely expanding ejecta. 

Subject headings: ISM: individual (Cassiopeia A) — supernova remnants — radiation mechanisms: 
non-thermal — acceleration of particles — shock waves 
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1. INTRODUCTION 

Supernova remnants (SNRs) are the main candidates 
for producing Galactic cosmic rays, with energies at least 
up to the so-called knee of the cosmic ray spectrum at 
~ 3x 10 15 eV. The first direct evidence for this is the de- 
tection of X-ray synchrotron emission ca used by ~ 10 14 
eV el ectrons (first established for SN1006. lKoyama et alJ 
[19951) . Since the energy of electrons suffers from radia- 
tion losses, this might indicate even higher energies for 
ions. Moreover, hard X-ray tails up to 8 keV have 
been discovered for several Galactic SNRs ([Alien et al.l 
1999). This has been contributed to either non-thermal 



bre msstrahlung (Laminell2001h or to synchrotron radia- 
tion (|Allen et alJ 1997t l . In recent years, additional direct 
evidence for efficient cosmic ray acceleration has come 
from detection of TeV 7-rays for several SNRs by the 
High Energy Gamma-Ray Astronomy (HEGRA) exper- 
ime nt, the High Energy Sp ectroscopic Syst em (H.E.S.S, 
e.g. lAharonian et al.lT2004h and MAGIC (jAlbert et all 
12007ft . The" 7-ray emission is either caused by inverse 
Compton scattering by the same electrons that cause X- 
ray synchrotron emission, or by pion production caused 
by collisions of accelerated ions with the background 
plasma. 

Cassiopeia A (Cas A ) is one of the sup ernova remnants 
with a hard X-ray tail (|The et al.lll996f) and has recently 
also been detecte d in 7-rays ([Aharonian et al.l 120011 : 
lAlbert et al.ll2007l ) . This remnant was until recently 1 the 
youngest known supernova remnant in the Galaxy; its 

1 Recent expansion measurements of Gl.9+0.3 show an 
age for this remnan t of around 100 years l(Grccn ct afj 120081 : 
Reynolds ct al.11200811 . 



super nova was probably around 1671 ([Thorstensen et ail 
2001). In 2001 Chandra detected thin, X -ray synchrotron 
emitting, filame nts at the forward sh ock of the remnant 
(Fig. [H see also iGotthelf et all 1200 fl ). This implies the 
presence of electrons with energies of ~ 10 13 eV for the 
magnetic fields in Cas A, estimated to be 0.1 mG t o 
0.6 mG (jVink fc Lamind[200l [Berezhko fe Volkll2004h . 
These synchrotron rims can be understood in the context 
of diffusive shock acceleration and synchrotron cooling 
downstream of the shock. 

In addition, the Chandra image shows thin filaments at 
the inside of the remnant. Some of these i nner filaments 
show a featureless spectrum (Fig. [5]). Hugh es et al.1 
(2000) identified one of the inner filaments at the West 
side of the remnant ('region D') as being the pro- 
jected forward s hock, based on its featureless spectrum. 
iDeLanev et al.l (|2004h found that the kinematics of the 
inner filaments, which they interpreted as projected for- 
ward shock filaments, are different from the forward 
shock; they have a lower velocity. 

Diffusive shock acceleration is a process which ac- 
cel erates cosmic rays at a shock front (for a review, 
see lMalkov fe DrurvH200il ). This mechanism accelerates 
charged particles of sufficient energy, which scatter on 
turbulent magnetic fields/plasma waves on both sides of 
the shock front. Each time the shock front is crossed, the 
particle gains energy, due to the difference in plasma ve- 
locity between both sides of the shock front. The higher 
the difference in the velocities is, the more energy is 
gained in one iteration and the higher the magnetic field 
and magnetic field turbulence, the more often particles 
cross the shock front. Since at the location where efficient 
particle acceleration takes place recently accelerated elec- 
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Fig. 1. — Million second 4-6 keV continuum image of Cas A, 
obtained with Cha ndra. The rectangle indicates region 'D' in 
IHuehes ct al. ( 2000T) . See the electronic edition of the Journal for 
a color version of this figure. 

trons are present, these locations show X-ray synchrotron 
radiation. However, further downstream from the shock 
front, synchrotron losses result in lower maximum ener- 
gies of the synchrotron radiation. 

The synchrotron spectrum can be approximated over 
a large range in frequencies with a power-law in flux 
density: F v oc v~ a with an index (a) related to the 
power-law index of the electron distribution (p) as: a = 
(p — l)/2. In what follows, we use index T, which 
refers to the number density index T = — (a + 1) and 
n(E) oc E v . Near the maximum electron energies, the 
electron spectrum has an exponential cut-off, but the 
resulting synchrotron spectrum cuts off less abruptly, 
rough ly as exp(— y/v/v max ) ( Zira kashvili fc Aharonianl 
2007). In contrast, the other important continuum emis- 
sion process, thermal bremsstrahlung, has an exponential 
cut-off (oc exp(— hv/kT)). In Cas A, t he plasma temper - 
ature ranges between 0.6 and 3.6 keV (|Yang et alj f2007). 
One should take into account that these temperatures 
originate from a thermal model. If a partly non-thermal 
spectrum is fitted with a thermal model, the fitted tem- 
perature tends to increase with respect to the real tem- 
perature of the plasma. For a thermal bremsstrahlung 
spectrum with a temperature of 3.5 keV, the power-law 
index between 4.2 and 6.0 keV is -2.8. For a synchrotron 
spectrum at the forward shock of Cas A, we typically 
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Fig. 2. — The spectrum of Cas A as observed by Chandra. 
Below is the spectru m of the featureless filament ('D') described by 
IHuehes ct al. (2000) extracted from the megasecond observation, 
above is a spectrum of the whole remnant of one single observation 
(ObsID 4638), multiplied by 0.1. 



measure a power-law index of -2.1. We therefore ex- 
pect the bremssstrahlung continuum to be steeper than 
the synchrotron continuum in the 4 to 6 keV continuum 
band. 

In this paper we investigate the shape of the continuum 
spectrum and its spatial distribution in order to address 
several questions pertaining to the shock acceleration in 
Cas A: What is the location of the X-ray synchrotron 
filaments? What fraction of the overall X-ray continuum 
is thermal and what fraction is non-thermal? And what 
are the implications for the hard X-ray emission, above 
10 keV. We do this by analyzing the Chandra megasec- 
ond observation of Cas A. 

2. THE USED DATA 

Chandra observed C as A for one million second in 2004 
from Februari to May (jHwang et alJl2004f ) with the ACIS 
S3. For extracting the images, we use the CIAO package, 
version 3.4. To align the separate pointings, the central 
compact object is taken as a reference point. We mostly 
concentrate on the 4.2 to 6 keV energy band, but in ad- 
dition, we also extracted an image in the line of Si XIII 
Hea from 1.80 to 1.90 keV. We corrected this for contin- 
uum emission by subtracting the average of two images 
next to the Si line: in the 1.63 - 1.64 and 2.13 - 2.14 keV 
bands. 

We compared the results of our analysis of the Chandra 
data with a VLA radio map, made in the 6 cm band 
in 2000-2001. Thi s map was kindly provided to us by 
Tracey DeLaney (|Delanev et al.l [20051 ) . The VLA image 
has a resolution of 0.38" x 0.33", which is comparable to 
the Chandra telescope resolution of 0.42". However, the 
Chandra pixels size of 0.49" slightly undersamples the 
Chandra resolution. 

In section HI we make use of BeppoSAX observations 
made in June 2001, with an exposure time of 501 ks. The 
hard X-ray spe ctrum obtained wit h the PDS instrument 
is described in IVink et all (|2001[ ) . The data obtained 
with the MECS instruments were never published before, 
but the deta ils of the analysis are similar to the analysis 
described in I Vink et al.l (|1999l) for a shorter exposure. 

Finally, the INTE GRAL-ISGRI spectru m used in sec- 
tion d] is described in lRenaud et al.l (|2006[ ). 

3. SEPARATING THE FORWARD AND REVERSE SHOCK 

The surface brightness (S(r)) of an optically thin ob- 
ject consists of the emissivity function of this object 
(e(r)), integrated along the line-of-sight. For a spheri- 
cally symmetric object, this integral is as follows: 

S(r) = 2 [ R e(r>) ^-^ dr 1 (1) 

In which R denotes the outer radius of this object. We 
do the deconvolution in cylindrical coordinates, with 8 
perpendicular to r, and 9 = is defined in the West, 
increasing counterclockwise. The number of counts in a 
sector at inner radius r, angular width d6 (in radians) 
and thickness dr is: 

C r = H{r)d0 (rdr + (2) 

By dividing Cas A in 18 sectors of 20° each, we first make 
surface brightness profiles of the X-ray continuum be- 
tween 4.2 and 6 keV. For the surface brightness profiles, 
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Fig. 3. — Radial surface brightness profile (smooth solid line) in 
the 4.2 to 6.0 keV energy band at an angl e of 10° t o 30° (includ- 
ing the featureless filament found by Hughes et al. (2000)). The 
smooth solid line indicates what the profile will look like if the sur- 
face brightness is produced by an emissivity function with only a 
peak at the outer shock. For the dashed line, we use an emissiv- 
ity function with two peaks. Note that this is not a fit, just an 
illustrative example. 



we adopt the center of expansion (|Thorstensen et ail 
[2001b : a = 23 h 23 m 27 s .77 and 5 = 58°48'49'.'4 (Equinox 
J2000) and take step-sizes of 0'.'5. Furthermore, we as- 
sume sperical symmetry for each sector individually. For 
an example of a surface brightness profile, see the ra- 
dial surface brightness profile between an angle of 10° 
to 30° in Fig. [3] In this figure, we can clearly see the 
outer shock coming up at 160", as already found by 
iGotthelf etaTI (|2001f ). Using equations [T] and H we now 
make an emission profile in such a way that it fits the sur- 
face brightness of the outer shock (smooth solid line) . As 
the line indicates, this thin, outer shell, can not account 
for all the surface brightness in the center of the remnant 
(the line which fits the outer peak is for a small R at least 
15 times lower than the surface brightness profile). We 
also see a peak in the surface brightness at 126". We can 
identify this surface brightness peak with the featureless 
filament described by Hughes et al. ( 2000). If we now in- 
clude a second emission peak at 126" (dashed line), we see 
that a large part of the surface brightness in the center 
is covered. 

To go from surface brightness to emissiv ity we use a 
general de-convolution method describe d bvlLucvi (119741) 
and p reviously used on SN1006 by IWillingale et al.l 
(1996). To test this algorithm, we simulated emis- 
sion functions, convolved them into a surface brightness 
profile and add Poisonian noise using the IDL routine 
' poidev' from the NASA IDL Astronomy User's Library 
(Landsman 1993), we made them in such a way that the 
peak of the number of counts in one bin of 20° x 0"5 
is 11000, which corresponds to the peak of Fig. [3] in 
terms of counts per 0'.'5 x 20° bin. After constructing 
this surface brightness profile we used the Lucy algorithm 
to recover the emissivity function. We stopped the de- 
convolution after 30 iterations. The resulting emissivity 
function is similar to the original one, within 10% for 
R > 40"and even within 5% for R > 75". 

3.1. Results of the deconvolution 

The results of the de-convolution of the Chandra data 
of Cas A are given in Fig. 2] We see the forward shock 
coming up in the X-ray continuum at 160"and with a 



width of 10". What we also see, is another ring inside 
the forward shock with a radius of 115"with a width of 
30", shifted to the right by 15". The forward shock and 
this inner ring are two almost complete circles, except 
for the NE region, where the jet is present. If we do 
the same for the radio and Si images of Cas A, we see a 
ring at the same location as the inner ring of the X-ray 
continuum. These inner rings in radio and Si have in 
previous resear ches been identified as the reverse shock 
(|Gotthelf et~aII[200lT) . 

We note that the results of the deconvolution are far 
from perfect; the xled °f the individual fits range from 
20 to 253 for 4.2 to 6.0 keV and from 18 up to 1150 for 
the Si band. This is to be expected, since the filamentary 
structure is clearly inconsistent with the assumed spher- 
ical sy mmetry. Moreover, th e X-ray Doppler maps of 
Cas A (jWillingale et al.ll2002f ) also show deviations from 
spherical symmetry. However, the fact that we can trace 
the inner and outer shells consistently from 18 indepen- 
dent deprojected sectors, argues for the veracity of the 
two rings. 

In order to estimate the total contribution of each shell 
to the total emissivity of Cas A in each band, we multi- 
ply the deprojected shells with the volume of the shells 
(4ttR 2 AR), taking into account the relative contributions 
of each sector. For the continuum band in the East to 
South (0 £ [180°, 260°]), for which the overall X-ray ra- 
diation is dominated by silicon line emission, 16 % comes 
from the outer ring and 84 % comes from the inner ring. 
For the Western part of the remnant (6 € [-100°, 20°]), 
in which most of the filaments dominated by continuum 
are situated, we find that 18 % of the continuum emis- 
sion is due to the outer ring and 82% due to the inner 
ring. 

4. THERMAL BREMSSTRAHLUNG VERSUS 
SYNCHROTRON RADIATION 

The hard X-ray e mission from Ca s A is clearly non- 
thermal in natur e (|The et all Il996t lAllen et all 119971 ; 
iFavata e~ al. 1997). In order to see how this relates to the 
Chandra spectrum of Cas A, we show in Fig. [5] both the 
broadband (1-8 keV) Chandra spectrum and the hard 
X-ray spectru m (above 20 keV ) as obtained with Bep- 
poS AX/PDS dVink et al.l [20ol|) and INTEGRAL/ISGRI 
(jRenaud et al.l 120061 ) . The model shown is fitted to the 
hard X-ray data only and is similar to the one used by 
Renaud et al., which includes the contributions of 44 Ti 
decay at 67.8 and 78.4 keV. We included absorption with 
Nh = 1-3 x 10 22 cm~ 2 . Of interest here are the best fit 
parameters of the non-thermal component, fitted with a 
power-law. Our best fit parameters are T = —3.4 ± 0.2 
and the normalisation is 3.22 ±1.9 counts keV -1 cm -2 
s _1 at 1 keV. We checked these results with previous 
res ults on the hard X-ray ta i l of t he spectrum of Cas 
A. iRothschild fc Lingenfelteil (|2003| ) observed Cas A for 
226 ks with RXTE and fitted a power-law to the HEXTE 
data from 20 to 200 keV. They find a power-law index of 
-3.125 ± 0.050 with a normalization of 2.0 ± 0.6 counts 
cm _2 s _1 at 1 keV. 

We see that the BeppoSAX power-law, extrapolated 
back to Chandra energies, fits the continuum between 
2.0 - 7.0 keV well. Within a confidence level of 90%, 
the extrapolated powerlaw contains a minimum of 93% 
of the counts in the Chandra continuum bands. The 
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Fig. 4. — In the left figure we see the de-convolution of the 4.2 to 6.0 keV continuum radiation. For all the figures, the deconvolution is 
done in sectors, within each sector, spherical symmetry is assumed. The scale is in % per bin. Each sector is scaled in such a way that the 
total adds up to 100%. The middle figure shows the de-convoluted Si image and the right figure the deconvolution of the radio emission. 
The emissivity is multiplied with the volume of the shell to get the total emission over the whole sphere for each dr. 
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Fig. 5. — Unfolded broadband X-ray spectrum of Cas A. The 
solid line represents the best fit for the combined BeppoSAX/PDS 
and INTEGRAL/ISGRI data. Extrapolation of this line back to 
the Chandra regime fits the continuum surprisingly well. The dot- 
ted line indicates the minimal number of counts in the 4.2 to 6.0 
keV band, allowed within the errors of the extrapolation of the 
hard X-ray fit. The Chandra broadband spectrum of Cas A was 
made from o ne of the observatio ns (ObsID 4638) of the megasecond 
observation (Hwang ct al. 2004). 

results of iRothschild fc Lingenfelterl ( 2003) amounts to 
a continuum flux from 4 to 6 keV which is 28% of the 
continuum flux measured by Chandra. Note that if we 
trust the power-law continuum model, there is little room 
for an additional thermal component. This is surpris- 
ing, since we know that the continuum in the 1-10 keV 
range is at least partially due to thermal bremsstrahlung, 
which must accompany the copious X-ray line emission. 
In fact, fitting the spectrum in the 0.5 -10 keV band can 
be done with a pure thermal model (|Vink et al.l 119961 : 
iWillingale et^MM) . 

To pinpoint the different mechanisms which contribute 
to the continuum in the Chandra band, we fit a power- 
law to regions of 4.9" x 4.9" (10 x 10 pixels) for the band 
between 4.2 and 6.0 keV. We use the eventfiles of the 
megasecond observation, which we merge, using the cen- 
tral compact object as a reference point. We use one 
ARF 2 file, made for the spectrum of Fig. [5l We es- 
timated the background contribution, using an annulus 
around Cas A, with R m i n = 205'.'8 and R m ax — 235'.'2 
and a center as defined in section [3l We neglected the 
effects of differential absorption over this small band; 

2 Ancillary Response File, which contains a description of the 
effective area of the instrument. 



in the most extreme cases, the absorption varies from 
1 x 10 22 to 1.7 x 10 22 cm -2 which has an effect of at max- 
imum 0.09 in the power-law index. Fitting the contin- 
uum of the total remnant between 4.2 to 6.0 keV, using 
Njj = 1.3 x 10 22 cm -2 , we found that the average index 
is -3.14. 

Fig. [6] shows that the power-law index varies consider- 
ably inside Cas A. We see that some of the regions with 
hard spectra overlap with known regions of non-thermal 
X-ray emission. For example the forward shock region 
appears to have a harder power-law. And also the non- 
thermal filaments in the West show up in the spectral 
index map. However, there are also is also hard contin- 
uum emission in regions where there is also line emission. 
This is illustrated by Figures [H] and [7] 

The spectral map of Fig [6] implies that if we observe 
Cas A at higher energies, the Western part of the rem- 
nant will become more prominent. We validated this 
by comparing an extrapolation of our best fit power-law 
spectra with the 501 ks BeppoSAX/MECS observation 
in the 9-11 keV band. 3 We deconvolved the BeppoSAX 
image using the Lucy deco nvolution method , following 
the procedure described in IVink et alJ (jl999) . The re- 
sults are shown in Fig [8] for comparison, we also show 
the 4 to 6 keV BeppoSAX image next to the Chandra 
image smoothed to roughly the same resolution (the ex- 
trapolation without smoothing is shown in Fig[9]). Quali- 
tatively the 9-11 keV image of BeppoSAX agrees with the 
extrapolated Chandra image: the South-Eastern part of 
the remnant is relatively less bright, whereas the Western 
part and the Southern part of the center are becoming 
more prominent. 

Since for a bremsstrahlung continuum, we expect an 
exponential cut-off and thus a soft spectrum and for 
synchrotron radiation, we expect a hard spectrum, we 
tentatively identify the hard spectra with synchrotron 
emission. Although there is a likely overlap in spec- 
tral index between thermal and non-thermal emission we 
can nevertheless estimate the total contribution of non- 
thermal emission by noting that in Fig. [6] left, those 
regions that have abundant line emission (the Eastern 
part of the shell), have power-law spectra steeper than 
-3.1. The total flux associated with those power-law in- 

3 There is also an observation made by XMM-Newton of Cas A 
up to 15 keV ( Blocker et al. 2001), because of errors in the exposure 
map of this observation, we choose to use the BeppoSAX data. 
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Fig. 6. — In the left Figure, we see the fitted power-law indices (r) + 3.1, to show the difference between the individual fits and the overall 
spectrum. Dark indicates a steeper spectrum. In the right Figure, we see a map of an image in the 4.2-6.0 keV continuum bands divided 
by a broadband image. The lighter color means re latively more continuum and thus a harder spectrum. Note the similarities between the 
images: as already noted by DcLancy et al. (2004), the continuum dominated areas tend to have a harder spectrum. Even so, not all hard 
spectra have a lack of line emission. 
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Fig. 7. — On the x-axis is the ratio of continuum over broadband 
in arbitrary units (the same units as in Fig. \E\ right image), a 
higher value means relatively more continuum. On the y-axis the 
fitted powerlaw indices of the corresponding bins. 

dices amounts to 46 % of the total flux in the 4.2 to 
6 keV band (Fig. HQ]). Suggesting that the other 54% 
of the flux, is due to non-thermal radiation. This cor- 
responds roughly to a total non-thermal flux above 4 
keV of 2.7 xl0 _10 erg s _1 cm" 2 . If we take a power-law 
index of -2.8 as an upper limit, this corresponds to a 
thermal bremsstrahlung model with kT = 3.5 keV, the 
non-thermal contribution comes down to 33 %. This is 
indeed a large fraction of the total contintuum emission, 
but not as large as the 93% suggested by extrapolation 
of the hard X-ray power-law. The reason for this is prob- 
ably that the true spectral shapes are not exactly power- 
laws, but steepen at higher energies. So apparently, in 
the Chandra band the addition of soft thermal emission 
and hard non-thermal gives by coincidence almost the 
same power law index, as the steepening non-thermal 
spectra at high energies. The lack of an obvious spectral 
break is probably due to the variation in power law in- 
dices for the non-thermal spectra, and the variation in 
cut-off energies across the remnant. 

4.1. West side of Cas A 

So far, we have shown that a large part of the non- 
thermal emission is associated with the reverse shock re- 
gion and that most of this emission is coming the West- 
ern part of the remnant. This prompted us to investigate 



Fig. 8. — The upper row is a Chandra image of Cas A in the 4 to 
6 keV band, convolved to the BeppoSAX resolution (left) and the 
BeppoSAX image of Cas A in the same energy band (right). The 
second row has on the left side the image of Fig. \E\ extrapolated to 
9 to 11 keV (Fig. |9j and convolved to the BeppoSAX resolution. 
The right image shows the actual BeppoSAX image in the 9 to 11 
keV band. The extrapolation matches qualitatively the BeppoSAX 
image and therefore, the powerlaw fits, give a reasonable prediction 
of the hard X-ray image of Cas A. 



the spectrum of the Western part more closely. Specifi- 
cally we were interested in whether the reverse shock in 
the Western part shocks the less dense material. The 
reason is that the Western region not only shows more 
non-thermal emission, but there also seems to be less 
thermal emission. A low density of a shocked material 
results, apart from less flux from shock-heated plasma, 
also in a low ionization time scale, n c t. With this in 
mind, we extracted spectra from one of the featureless 
filaments, located at the reverse shock. T his is the same 
filament as used by iHughes et alj f2000) (filament 'D'). 
We extracted the spectrum of this region in the individ- 
ual eventfiles of the megasecond image and added them 
with the addspec tool of ftools in heasoft. We fit the 
resulting spectrum with a model consisting of a power- 
law and the non-equilibrium ionization (NEI) model of 
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Fig. 9. — An extapolation of the best fit power-laws from Fig. [6] 
to the 9 to 11 keV band 
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Fig. 10. — A plot of the cumulative flux between 4.2 and 6 keV 
versus the best-fit powerlaw index as found in Fig [6] The solid 
line shows that 46% of the flux has a power-law index steeper than 
-3.1. The dotted line shows that 33% of the flux has a power-law 
flatter than -2.8; the power-law index associated with a thermal 
bremsstrahlung model with kT = 3.5 keV. 



TABLE 1 

Results of the spectral fit to the featureless 
filament 



Parameter 



Value 



n H (cm 2 ) 
n e t (cm -3 s) 
kT (keV) 

r 

7 P ow/F nci (E 6 [0.4 keV,8 keV]) 
X 2 /clc,f 



10 1 



-0.03 

3±i:8 >< 

2+0.15 
'-0.15 
(-0.07 
-0.07 
2.23 
283/145 



I.. 
-2.09^ 



Note. - Prev ious studies ( Kcohanc 119961 : 

Willingalc et al. 2002) show that the absorption in the 
West side is higher than in the East, with values up 
to 1.7 xlO 22 cm -2 . So a value of 1.67 xlO 22 cm~ 2 is 
still within the expectations. 

SPEX (|Kaastra et al.|[l996f ). The best fit parameters are 
listed in Tabled 

5. DISCUSSION 

Using the a combination of deprojection and character- 
ization of the X-ray continuum power-law slope of the one 
megasecond Chandra observation, we have determined 
the location of the non-thermal X-ray emission of Cas A. 
We have found that a significant fraction of the X-ray 
continuum emission is coming from two separate rings; 



one can be identified with the forward shock, the other 
with the shell dominated by bright ejecta. About 54% of 
all continuum emission in the 4.2 to 6 keV band is likely 
to be non-thermal. 

The non-thermal emission from the forward shock is 
likely due to synchrotron radiation from electrons ac- 
celerated at the forward shock, confine d to the forward 
shock region by synchrotro n cooling (|Vink fc Lamina! 
I2003t iBerezhko fc Volkj|2004h . However, the non-thermal 
emission of Cas A as a whole is dominated by the contri- 
bution of the inner ring. This inner region of emission is 
shifted to the West of the remnant by ~ 15". A similar 
shift was found for the location of the r everse shock by 
iReed et all (|1995ft ; lGotthelf et alj (|2001h . on the basis of 
optical data and a deprojection of radio data and a 50 
ks Ch andra observation, in the Si band. iGotthelf et al.l 
( 2001) also found marginal evidence for a higher emis- 
sivity at the reverse shock in the 4 to 6 keV band using 
the 50 ks Chandra observation. 

5.1. Nature and location of the inner non-thermal 
emission 

Since some of the spectra in the West hardly show 
any line emission, we think it is likely that this is 
synchrotron emission, since non-thermal bremsstrahlung 
without line emission is hard to establish, unless one 
has peculiar abundances. Such a bremsstrahlung mode l 
was once invoked for SN1006 (jHamilton et al.l fl986). 
but has now been abandoned in favor of synchrotron 
emission. Furthermore, a recent paper by IVmkl (2008) 
shows that electrons with energies close to the ther- 
mal electron energy distribution, lose their energy rel- 
atively fast due to Coulomb losses. This process is al- 
ready important in the Chandra energy band for n c t < 
8.6 x 10 10 cm _3 s; the value reported here for a filament, 
dominated by continuum emission in the West (section 
I4TT] Table [Q. For n e t ~ 10 n cm~ 3 s, typical for Cas 
A (jWillingale et al.l [20021 ). one only expects to see non- 
thermal bremsstrahlung for photon energies > 100 keV. 

When it comes to the location of the synchrotron emis- 
sion, there are two possibilities: the reverse shock and the 
contact discontinuity. The contact discontinuity marks 
the border between shocked ejecta and shocked circum- 
stellar medium. Hydrodynamical solutions show that 
the density and the magnetic field peak at this radius 
(|Chevalierl Il982t ILvutikov fc Pohll I2004D . For a super- 
nova remnant evolving in a stellar wind, the contact dis- 
continuity is close to the reverse shock. So from our 
deprojections it difficult to judge whether the X-ray syn- 
chrotron emission is coming from the reverse shock or the 
contact discontinuity. It is unlikely that electrons are ac- 
celerated to TeV energies at the contact discontinuity, 
since no viable acceleration me chanism is known (how- 
ever, see ILvutikov fc Po hl 2004), but there are two ways 
of generating X-ray synchrotron emission at the contact 
discontinuity: 1) Due to an increase of the magnetic field, 
electrons with relatively low energies suddenly light up in 
X-rays. 2) High energy electrons and positrons are cre- 
ated through the decay of charg ed pions (n^) , caused by 
hadronic cosmic ray collisions (G aisser l 1990). Charged 
pions decay into muons and muon neutrinos. The muons 
decay into electrons/positrons and electron and muon 
neutrinos. The electrons and positrons thus created are 
often called secondary electrons and positrons. 
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For option 1: For synchrotron radiation the relation 
between electron energy and photon energy is 

E ph = 19E* eV B G keV (3) 

(jGinzburg fc Svrovatskiilll965f ). For electrons to be in- 
visible in X-ray synchrotron radiation, the peak energy 
should be an order of magnitude lower than the X-ray 
continuum at 4-6 keV. Therefore, the magnetic field 
should increase an order of magnitude at the contact 
discontinuity. As an example, an electron, accelerated at 
the forward shock, with an energy of 3 TeV in the typ- 
ical magnetic field of Cas A of 0.5 mG, typically emits 
photons of 0.09 keV. When this electron suddenly en- 
ters a region of 5 mG, it will emit photons of typically 
1 keV; detectable in X-rays. However, the synchrotron 
loss time for such an electron in a 0.5 mG magnetic field, 
18 years, is rather short compared to the age of Cas A. 
In this time the particle has to diffuse from the forward 
shock region to the contact discontinuity, wh ich is about 
0.5 p c, for a distance to Cas A of 3.4 kpc (Ree d et all 
I1995H . For 3 TeV and B=0.5 mG the diffusion constant 
is D = i] ■ 2 x 10 23 cm 2 s -1 , with n = 1 corresponding to 
Bohm diffusion (|Malkov fc Drurvll200lD . In 18 yr a par- 
ticle can diffuse by R w \/2Dt = 5 x 10~ 3 ^/r] pc. In other 
words, this model only works if the magnetic field turbu- 
lence is very low, corresponding to 77 > 10000. This un- 
likely, since the ample presence of cosmic rays results in 
magne tic field turbulence. Moreove r, at the shock front 
77 ~ 1 (|Vinkll2006l : rStage et al.ll2006h . For low values of 77, 
the diffusion length (Zdiff the length for which advection 
dominates over diffusion), is short: Z^iff = x 10~ 4 
parsec. So, for low 77, the advection velocity, (u = V[ s /4, 
with Vf s the forward shock velocity) is the relevant ve- 
locity. For a Vf s of 5800 kins" 1 , u — 1450 kms" 1 . 
Using this velocity, it takes 674 years to go from the for- 
ward shock to the contact discontinuity, longer than the 
synchrotron loss time of 18 years. 

It has been argued that the magnetic field at the 
shock front is hi gh, b ut rapidly decays towards the in- 
side dPohl et a l. 2005), increasing the synchrotron loss 
times, and increasing the diffusion constant. However, 
the decay in magnetic field should be reflected in the ra- 
dio emissivity from the forward shock to the inside, which 
is contrary to observations, that show gradual increase in 
emissivity toward th e center starting at the shock front 
dGotthelf et al.ll200l . 

For option 2 (secondary electrons): if the X-ray syn- 
chrotron emission from the inside is due to secondary 
electrons, this would be an important discovery. It would 
be evidence for the presence of TeV ion cosmic rays. 
Neutral pions (w ) are made in comparable quantities to 
charged pions. The power in secondary electrons should 
therefore be comparable to the luminosity in pion decay. 
In order to compare the fluxes from X-ray synchrotron 
and 7-ray radiation one should take into account the 
conversion from pion energy (it & 7T ) to photon en- 
ergies, both in the TeV band (through tt° decay) and 
in the X-ray band, due to synchrotron radiation from 
secondary electrons. Taking into account these various 
decay channels, one finds that E p h ~ 20£ 2 ToV _B G keV 

4 , so the X-ray synchrotron flux above 4 keV should 
correspond to roughly the 7-ray flux above 20 TeV. 

4 The decay product of n° decay is two photons, with each, 



We found for the X-ray synchrotron flux above 4 keV 
2.7x 10~ 10 erg s^'cm" 2 . This is a factor 130 higher than 
the 7-ray flux above 1 TeV, which we calculated to be 
2.1 x 10 -12 erg cm" 2 s _1 using the photon-flux and pho- 
ton in dex above 1 TeV re ported by the MAGIC collabo- 
ration (I Albert et al.l2007fl . Since we should actually have 
evaluated the flux above 20 TeV, we can rule out that 
the X-ray synchrotron emission is caused by secondary 
electrons. 

In our view it is therefore most likely that the X-ray 
synchrotron emission is caused by electrons accelerated 
at the reverse shock. Similarly to the forward shock re- 
gion the electrons are likely to be confined to a region 
near the shock itself. For a long time the reverse shock 
as location for acceleration has been neglected, because 
of its putative low magnetic field. Moreover, the abun- 
dance pattern of cosmic rays is consistent w ith accel- 
eratio n from plasmas with solar abundances (|Horandell 
2008). As far as the ma gnetic field is concerned, how- 
ever, lEllison et alj (2005) h ave argued that if magnetic 
field amplification (|Bell fc L ucck 2001) works at the for- 
ward shock it is likely to operate at the reverse shock as 
well. Also for the SNR RCW 86 it has been suggested 
that the X-ray synchrotron emission in the Southwest of 
the remnant i s coming from el ectrons accelerated at the 
reverse shock (jRho et alj |2002). 

5.2. Reverse shock velocity in the West and the 
presence of synchrotron emission 

The presence of X-ray synchrotron radiation from 
shock accelerated electrons is only expected if the 
shock velocity is high enough: shock acceleration 
theory predicts for the maximum photon energies 
(|Aharonian fc Atovanlll999D : 



E, 



P h 



O.577" 



• 2000 kms" 



kcV, 



(4) 



with v s the shock velocity. Note that the photon en- 
ergy is independent of the magnetic field. For the re- 
verse shock the velocity in Eq. prefers to the shock speed 
in the frame of the ejecta. The ejecta velocity is equal 
to the free expansion velocity Wf. j — r c -Jt. For the re- 
verse shock the shock velocity as seen by the ejecta is 
therefore v s>e j = Uf, C j — w s ,obs, with f Si0 b s the shock ve- 
locity in the frame of the observer. The presence of 
X-ray synchrotron radiation from the reverse shock in 
the Western half of Cas A suggests a higher reverse 
shock than in the rest of the remnant. Proper mo- 
tions of knots at the insid e of Cas A in X-rays were 
most recently measured by iDeLanev et al.1 (|2004h . but 
some details of the measurements, including measure- 
ments of the proper motions as a function of azim uth, 
only appeared in appendix 4.4 of iDelanevI (|2004 Fig. 
4.6). In this Figure, we see that in the West the ex- 
pansion rate is between -0.1%yr _1 and Q.1%" 1 this cor- 
responds to i> r , b at the reverse shock of approximately 

on average, an energy E 7 ToV = 0.5E n o- For the charged pions, 
7T , the final decay product consists of electrons (positrons) and 



neutrinos ( 7r d 



+ fe(fe) + V^Vfj)). In 



this case the final electron or positron also take up, on average, 
half the initial pion energy E e ± fs 0.5E n ± ; hence, for a given pion 
energy, E e ± f» E 7 TeV ■ This electron or positron emits synchrotron 
radiation at a peak frequency given by Eq. \3\ Hence the close 
resemblance of this equation with Eq. |5| 
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-2000 km s 1 to 2000 km s 1 , implying a shock velocity 
5 u s , e j ~ 3900 - 79 00 kms~| In t he Northern and East- 
ern part of Cas A, iDelanevi (|2004h finds expansion rates 
of 0.2%yr _1 , corresponding to v t ,ob ~ 4000 kms" 1 and 
a shock velocity of u s , C j ~ 1900 kms" 1 at the reverse 
shock. 

In the optical iMorse et alj (|2004f ) reports a shock ve- 
locity of v s .cj ~ 3000 kins" 1 for a filament in the North- 
west. In our maps (Fig. [S]left) we find that this region 
has a hard power law slope, suggesting the presence of 
synchrotron radiation, but the overall X-ray emission is 
dominated by thermal X-ray line emission (Fig. [6] right). 

It indeed looks like the presence of X-ray synchrotron 
emission from the inner ring corresponds with reverse 
shock velocities in excess of 2000 kms" 1 . To some ex- 
tend this is surprising, since the X-ray synchrotron cmis- 
sivity function is rather broad and can result in some syn- 
chrotron emission beyond the maximum photon energy 
as defined in Eq.[H On the other hand, rj = 1 corresponds 
to Bohm diffusion, and represents the case for maximum 
acceleration efficiency. The fact that 2000 kms -1 seems 
close to the velocity dividing the presence or absence of 
X-ray synchrotron radiation again suggests that acceler- 
ation takes place close to the Bohm limit. Equation |4] 
suggests that the spectra in the West of the remnant will 
be harder than in the rest, where i> SiC j is lower. 

Note that the high reverse shock velocities in the West 
(and low observed velocities) are not in agreement with 
analytic hydro dynamic solutions for an SNR evolving in 
a stellar wind ( Laming k. H wang 2003). Adapting the 
parameters of this model such that the forward shock 
radius and velocity and the reverse shock radius, match 
those of Cas A, we find w SjC j ranging between 1600-2300 
kms" 1 , significantly lower than in the West, but in 
agreement with those in the Eastern part of Cas A. This 
suggests that the circumstellar density structure is more 
complex in the West. 

6. CONCLUSIONS 

We have presented an anaysis of the spatial and spec- 
tral variation of the X-ray continuum emission of Cas 
A, based on the 1 Ms Chandra observation. We find 
that harder continuum spectra are associated with the 
filaments, dominated by continuum emission, suggest- 

5 We have checked this result using the Chandra 1 Ms observa- 
tion and an obser vation of Chand ra made in 2000 using the proce- 
dure described in Vink ct alj 1119981 ). We concentrated on filament 
'D'. The result confirms the lower expansion or even backward ve- 
locities: we found —970 ± 140 kms" 1 , i.e. the filament seems 



ing that the harder spectra are caused by non-thermal 
radiation. A dominant fraction of the non-thermal 
emission appears to come from the reverse shock re- 
gion. We have discussed various options for the na- 
ture of the non-thermal emission and its origins. Some 
of our conclusions were indep endently also obtained by 
lUchivama &: Aharonianl (|2008[ ) , but based on X-ray vari- 
ability. Based on our analysis and discussion we come to 
the following conclusions: 

• The power-law index of the spectrum between 4.2- 
6.0 keV is an indicator for X-ray synchrotron emis- 
sion: there is a correlation between filaments, dom- 
inated by continuum emission and hard spectra, 

• hard X-ray spectra are not exclusively associated 
with filaments, dominated by continuum emission, 
suggesting that non-thermal emission comes also 
from other regions, 

• the non-thermal X-ray emission is likely to be syn- 
chrotron radiation, 

• the non-thermal accounts for about 54% of the 
overall continuum emission in the 4-6 keV band, 

• in the Western part of Cas A, most X-ray syn- 
chrotron comes from the reverse shock, 

• the dominance of X-ray synchrotron emission 
from the West is probably the result of a lo- 
cally higher reverse shock velocity of v s ~ 
6000 kms -1 (corresponding to a lower proper 
motion) than in the Eastern region (v s ~ 
1900 kms" 1 ). 
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to move towards the center. This implies a v s 



6900 kms" 



IVink et alj i 19981) already reported a lower expansion for the whole 
Western part of Cas A. 
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